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We present a class of models that describe self-diffusion on 
FCC(OOl) metal substrates within a common framework. The 
models are tested for Cu(OOl), Ag(OOl), Au(OOl), Ni(OOl) and 
Pd(OOl), and found to apply well for all of them. For each of 
these metals the models can be used to estimate the activation 
energy of any diffusion process using a few basic parameters 
which may be obtained from experiments, ab-initio or semi- 
empirical calculations. To demonstrate the approach, the pa- 
rameters of the models are optimized to describe self-diffusion 
on the (001) surface, by comparing the energy barriers to a full 
set of barriers obtained from semi-empirical potentials via the 
embedded atom method (EAM). It is found that these mod- 
els with at most four parameters, provide a good description 
of the full landscape of hopping energy barriers on FCC(OOl) 
surfaces. The main features of the diffusion processes revealed 
by EAM calculations are quantitatively reproducible by the 
models. 



I. INTRODUCTION 

Thin film growth processes involve complicated kinet- 
ics giving rise to a rich variety of surface morphologies. 
Within this vast domain, the study of the growth in 
the submonolayer regime is of particular interest due to 
the large impact of the initial kinetics on the resulting 
film structure. Experiments on thin film growth on well 
characterized substrates using molecular-beam epitaxy 
(MBE) have provided a large body of information about 
growth kinetics and morphology, and revealed that for 
a variety of systems and a broad temperature range, is- 
land nucleation is the dominant mechanism for crystal 
growth ||l|,0[. Diffraction methods such as helium beam 
scattering 10,^-^ , low energy electron diffraction 
and other techniques jri| , p^ , provide information on the 
collective behavior and the statistical properties of the 
surface. These techniques have been used to measure the 
island size distribution, the island density, and their scal- 
ing properties with respect to the coverage and the flux 
1^,0 . The variation of the island density with respect 
to the temperature was also studied 

More detailed information at the atomic scale is pro- 
vided by scanning tunneling microscopy (STM) P,p^-^ . 
Most notably, STM provides means to study the vari- 
ety of morphologies encountered in the different systems, 
or in the same system under different growth conditions 



[Q,|2-20 2^. In some experiments STM was used to ac- 
quire information on larger scales, e.g., island size dis- 
tributions [|l2|,^ 20 - 2^,Q . Despite the wealth of exper- 
imental results at the atomic scale, for decay rates of 
small islands, mobility of small islands and edge diffu- 
sion |^ , p3| - p7| , p9pl| ] , the underlying energetics is mostly 
inaccessible to direct experimental measurements. Thus, 
one must rely on theory to extract activation energies 
from the experimental results, and these are usually lim- 
ited in number, and sometimes are subject to alternative 
interpretations. 

The only technique which provides direct access to dif- 
fusion processes and activation energies at the atomic 
scale is field ion microscopy (FIM) |3^-|3q]. This tech- 
nique was used to identify the diffusion modes of adatoms 
|2|,|| as well as small islands [|4|j36) on FCC (001) metal 
surfaces, to measure their diffusion coefficiGnts, and to 
determine the sticking process of adatoms to an island 
psf . Recently there were several attempts to use STM 
to derive such local information directly [ p4| . 

Theoretical studies aimed at providing better under- 
standing of the relation between key processes at the 
atomic scale and the resulting morphologies have been 
done using Monte Carlo (MC) simulations [^-^. In 
simulations of island growth during deposition, atoms 
are deposited randomly on the substrate at rate F [given 
in monolayers (ML) per second] and then hop, attach 
to, and detach from existing islands according to some 
model. A common approach is to assume some key pro- 
cesses and their rates, and then simulate the growth pro- 
cess 38 4^,Q . In some cases information such as dif- 
fusion length, typical distance and time between nucle- 
ation events is assumed to be known a-priori, and is put 
by hand into the simulation in order to accelerate the 
computation . The advantage of this approach 

is that the models are well defined and use only few pa- 
rameters. These models are useful for studies of scal- 
ing and morphology but cannot provide a quantitative 
description of diffusion on a particular substrate. Fur- 
thermore, they account only for a limited number of pro- 
cesses, that are assumed to be the only significant ones. 

A complementary scheme employs the underlying ac- 
tivation energies. In this scheme the hopping rate h (in 
units of hops per second) of a given atom to each unoc- 
cupied nearest neighbor (NN) site is given by 

h = iy cxp{-EB/kBT) (1) 

where ly = 10^^ s^^ is the commonly used attempt rate. 
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Eb is the activation energy barrier, ks is the Boltzmann 
constant and T is the temperature. 

The activation energy barrier Eb depends on the local 
environment of the hopping atom, namely the configu- 
ration of occupied and unoccupied adjacent sites. Two 
approaches have been taken in the construction of the 
energy barriers for hopping in the simulations. One ap- 
proach was to construct simple models that include the 
desired features, such as stability and mobility of small 
islands, and that take into account properties such as 
bond energies |p|,p|,p7p9| |5^,p^,p6[ . In general, this ap- 
proach encompasses both the virtues and the drawbacks 
of the simpler approach presented before. 

A second approach is based on the use of an approxi- 
mate many-body energy functional to calculate the hop- 
ping energy barriers for a complete set of relevant con- 
figurations |^,^,|3 53 S^JI^ . This approach provides a 
good description of diffusion processes on the given sub- 
strate but only limited understanding due to the large 
number of parameters. 

In this paper we extend and further explore a frame- 
work for a systematic derivation of simple models for 
self-diffusion on FCC(OOl) surfaces out of a detailed and 
complicated set of energy barriers. Simple in this con- 
text means that only a small number of parameters are 
involved and all have a definite and intuitive interpreta- 
tion. Using sensible assumptions about the bond ener- 
gies and diffusion paths we obtain simple formulae for 
the activation energy barriers. We then optimize the pa- 
rameters of these formulae for each metal separately by 
using energy barriers obtained from the embedded-atom 
method. This procedure gives rise to simple models that 
have at most four parameters and provide good quan- 
titative description of the landscape of hopping energy 
barriers. In a previous publication we have introduced 
the framework and apphed it to Cu/Cu(001) growth [ p9| . 
Here we make a three-fold step forward, (a) We include 
four other FCC metals in the model and derive the ap- 
propriate parameters for them. This shows the utility of 
the models and provide a unifying framework that applies 
to a large class of metals; (b) Non-linear interactions are 
introduced in addition to the linear interactions consid- 
ered before. This allows for a more accurate optimization 
without increasing the number of parameters; (c) We ex- 
plore the basis of the assumptions underlying this scheme 
and the extent of their applicability. 

The paper is organized as follows. In Sec. II we intro- 
duce the physical framework of the model and discuss its 
underlying assumptions. The results of the EAM calcu- 
lations are given in Sec. III. The models are introduced 
in Sec. IV with the fitting to the EAM resuffs. This is 
followed by a discussion of the results and their implica- 
tions in Sec. V. 



II. APPLICABILITY CONSIDERATIONS 

The framework developed in this paper assumes several 
characteristics of the diffusion processes considered. It 
applies to systems for which these assumptions are valid, 
which includes most of the FCC metals in the moder- 
ate temperature regime (« 200 — 500 AT). The following 
assumptions are employed throughout the discussion. 
Bridge-site hopping of adatoms is in general dominant 
over exchange hopping (Fig. There has been a contro- 
versy concerning this assumption. Using semi-empirical 
methods the barrier for exchange hopping for Cu(OOl) 
was estimated to be 0.2 eV [po[-^2|, in agreement with 
the experimental data that was available at that time 
This is much lower than the barrier for bridge-site 
hopping. However, theoretical work using several other 
methods []63|-^ including EAM indicate that the ex- 
change barrier for Cu is much higher (more than 0.8 eV) 
and therefore bridge-site hopping is dominant. This con- 
clusion is also supported by recent experimental work and 
reinterpretation of the previous findings In general, 
bridge-site hopping is found to be dominant in all the 
metals studied here, (for Au, some exchange processes 
appear to be significant, yet in most cases bridge hopping 
is favorable). Due to the exponential dependence of the 
rate of each process on the corresponding energy barrier, 
it is generally reasonable to take into account only the 
mechanism which is energetically favorable (bridge-site 
hopping, in this case) and neglect the mechanism which 
exhibits higher activation energy barrier (exchange mech- 
anism) . 

Only nearest and next-nearest neighbor interactions are 
significant. Within this assumption one can obtain the 
activation energies for most diffusion processes to a good 
accuracy. However there are some processes such as va- 
cancy diffusion, where a larger environment affects the 
diffusing atom. 

There is one common attempt frequency for all processes. 
Since there is no systematic knowledge about the depen- 
dence of the attempt frequency on the local environment, 
the assumption of one common frequency is the usual 
practice. Estimations for the attempt frequency can be 
obtained using MD simulations |Q. Another way is to 
fit molecular static (MS) data to harmonic potential to 
find an effective force constant, and then deduce the fre- 
quency of oscillations. So far, little work has been done 
on this subject. Previous works, including interpretation 
of experimental results, usually pre-suppose some com- 
mon attempt frequency in the range 10^^-10^^s~^ []66[. 



III. THE EAM BARRIERS 

The models described in this work are tested by fit- 
ting their parameters to energy barriers for self-diffusion 
of Cu(OOl), Ag(OOl), Au(00]_K Ni(OOl) and Pd(OOl) sur- 
faces, obtained using EAM This method uses semi- 
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empirical potentials and provides a good description of 
self-diffusion on such surfaces Specifically, for all 

the metals considered here the EAM functions developed 
by Adams, Foiles, and Wolfer (AFW) ||6^ are employed. 
These functions are fitted to a similar but more accu- 
rate data base as the one employed by Foiles, Baskes, 
and Daw The calculations are done on a slab of 20 
square layers with 100 atoms in each layer. 

When an atom on the surface hops into a vacant near- 
est neighbor site it has to cross the energy barrier be- 
tween the initial and final sites. We have used molecular 
statics in conjunction with the EAM functions to find 
that energy barrier. This is simply the difference be- 
tween the energy at the bridge site (or more precisely, at 
the point along the path with highest energy) and in the 
initial site. 

The hopping energy barriers are calculated for all local 
environments as shown in Fig. |^, where seven adjacent 
sites, « = 0, . . . , 6 are taken into account, according to the 
assumptions presented in Sec. IL Each one of these sites 
can be either occupied {Si = 1) or vacant {Si = 0), giving 
rise to 2^ = 128 barriers. A binary representation is used 
to assign indices to these barriers. For each configuration 
{So, . . . ,Se) the barrier is given by Eg, where 

6 

n = Y,S^■2' (2) 

i=0 

takes the values n = 0, . . . , 127. The full set of hop- 
ping energy barriers (given in eV) is presented in Table 
I for Cu(OOl), Ag(OOl), Au(OOl), Ni(OOl) and Pd(OOl). 
To show these values in a compact form, each barrier in 
Table I corresponds to a configuration in which the oc- 
cupied sites are the union of the occupied sites in the 
picture on top of the given column and on the left hand 
side of the given row. The column in Table I in which 
a given configuration appears is determined by the occu- 
pancy of sites i = 2, 3, 6 while the row is determined by 
sites i = 0, 1, 4, 5. One can define 

1=2,3,6 1=0,1,4,5 

such that for each configuration n ~ ni + 712- To demon- 
strate the use of Table || , we will check for Cu(OOl) the 
barrier of the configuration in which sites 0, 3 and 4 are 
occupied and all other sites adjacent to the hopping atom 
are vacant. For this configuration, according to Eq. (||), 
ni = 8 and 712 = 17 {n — 25). The barrier, that is found 
in the column with the index 8 and the row with index 
17, in the line of Cu, is Ef = 0.89 eV. 

In Table ^ we use the symmetries of the configurations 
in the 3x3 cell (Fig. |^) to reduce the number of en- 
tries. There is a mirror symmetry plane perpendicular 
to the surface and containing the arrow of the hopping 
atom. Consequently, the columns of ni = 4 and 12, in 
which site j = 2 is occupied, stand also for the symmetric 
configurations in which i = 6 is occupied. In the other 



four columns, there are some configurations that, due to 
symmetry, appear twice. In such cases, the barrier for 
the configuration with larger n appears in italics. 

For the purpose of the calculations and parameteri- 
zation of the model we consider only hopping moves in 
which a single atom hops each time. It turns out, how- 
ever, that in some cases the molecular statics calcula- 
tions, used to obtain the barriers, give rise to concerted 
moves. In such moves the atom at site i = 3 follows the 
hopping atom and takes the place vacated by the hop- 
ping atom. This fact significantly reduces the barrier. 
It turns out that for configurations in which concerted 
moves appear, they can be suppressed by adding a col- 
umn of three atoms on the left hand side of sites i = 0, 3 
and 4. In Table I, the energy values for those configu- 
ration in which a concerted move was found, are shown 
in parenthesis. The barrier obtained when the concerted 
move was suppressed is shown to the left of the paren- 
thesis. 

To gain a better understanding of the barrier en- 
ergy landscape we present the barrier height distribu- 
tion (without concerted moves) in Fig. ^ for the five 
metals considered. We observe that this distribution ex- 
hibits four groups. This feature is in agreement with Ref. 
psf where a different method was used to calculate 
the barriers. Each group, corresponds to a single or a 
double column in Table I. In general, group I includes 
very fast moves towards island edges, group II includes 
moves along the edge, group III includes, most notably, 
the single atom move , while group IV includes detach- 
ment moves. 



IV. THE MODELS 
A. The Additivity Assumption 

The starting point in the construction of a simple 
model that describes the hopping energy barriers for all 
the configurations of Fig. ||, is the assumption that the 
contributions of all adjacent atoms to the energy barrier 
add up linearly. To examine this assumption for Cu(OOl), 
we evaluated directly the binding energies within the 
EAM approach for a series of configurations from which 
we extracted the relevant bond energies. 

The binding energy of a given configuration of adatoms 
on the surface is evaluated as follows: First we calcu- 
late the total energy of the system in that configuration. 
Then we find the total energy of another configuration in 
which there is the same number of adatoms on the sur- 
face but they are far apart from each other, (by that we 
mean that moving them one lattice site in any direction 
would not change the total energy). The binding energy 
between the adatoms is given by the difference in the to- 
tal energies between the two configurations. An example 
of the procedure is shown in Fig. |j. It appears, that the 
evaluation of the NNN bond energy is easier since one can 
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construct a sufRciently large set of configurations which 
include only NNN bonds with no NN bonds. In the case 
of NN bonds, most of the relevant configurations also 
include NNN bonds [Fig. || (a)]. Similarly, considering 
an atom on top of a bridge site, typical configurations 
which include atoms adjacent to the bridge site exhibit 
NN bonds between them as shown in Fig. ||(b). 

Therefore, we will first examine the additivity of the 
NNN bonds employing a series of four configurations in 
which an adatom has 1, 2, 3 and 4 NNN on the sur- 
face. For each one of these four configurations the total 
energy is compared to that of a configuration with the 
same number of adatoms, in which they are far apart 
from each other. In Fig. ^(a), the total binding energy 
between adatoms is plotted as a function of the number 
of NNN bonds. The best linear fit is drawn, and its slope 
yields a value of E^nn — 0.0512 eV. The next step is to 
examine the linearity of the NN binding energy using a 
series of four configurations in which an adatom has 1,2, 
3 and 4 NN's on the surface, within a similar procedure. 
The NNN bonds in each configuration are deducted using 
the value obtained before. The results are shown in Fig. 
^(b) and the NN bond energy is obtained: E^n = 0.324 
eV. A similar analysis for an adatom on a bridge site is 
shown in Fig. ^(c) and the binding energy between an 
atom on the bridge site and an adjacent atom is given by 



E 



NN{bridge) 



0.345 eV. 



B. Construction of the Models 

The energy barrier Eb for a certain process is the 
difference between the binding energies of the hopping 
adatom (to the substrate and to adjacent adatoms) at 
the initial position, i?i„, and at the bridge site, Etop, 
namely Eb = Etop — Em- On the basis of the additiv- 
ity feature just demonstrated, we will now express these 
binding energies as the sum of the occupation states of 
the relevant sites. The first approximation for the ener- 
gies gives a model (model I) with only two parameters, 
that reproduces the main features of the EAM barriers. 
In order to establish the model, there are t wo th ings to 
note about the parameters obtained in Sec. EVA. First, 



the values of NN binding energies at the lattice site and 
bridge site are very close. This reflects the fact that the 
NN distance corresponds approximately to the minimum 
potential of the two-body interaction. Second, both these 
energies are much larger than the NNN binding energy. 
These two features are quite general and common to all 
the metals we discuss here. For the simplest model we 
will neglect the effect of the NNN atoms, and assume a 
single NN binding energy, AE^n, for both lattice and 
bridge sites. The resulting expression for the binding en- 
ergy at the initial (fourfold hollow) site is 

Ell^ = El - AEnn -{S. + Ss + S^). (4) 

The energy of an isolated atom is E^ . The energy of the 
hopping atom when it is on the bridge site is given by: 



^top ~ ^top ■ 



AEj\iN • (5*1 -l- 5*2 + S5 + Sf) 



(5) 



where E^^p is the energy of an isolated atom on top of a 
bridge site. Thus, for a given configuration the barrier, 
E"^ — E'j'Jj, for an atom to hop into an adjacent 



vacant site is given in model I by: 



- 



AEnn ■ {S3 — S2 — Se) 



(6) 



where E°, = E^^ 



E^ and n is given by Eq. (j|). In this 
model only three sites affect the energy barrier, which can 
take only four different values, as the expression in the 
parenthesis can be either 1, 0, -1 or -2. Each of these 
four barrier values corresponds to one of the four groups 
in Fig. ||. The parameters of this model, as well as those 
of the models discussed below, are adjusted to best fit 
the EAM data. More specificly, we found the parameters 
that best describe the 128 EAM barriers by minimizing 
the sum of squares: 
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R = ^[E'^iEAM] - E'^iModel)]^ 



(7) 



n=0 



The values obtained for these parameters for the five met- 
als are shown in Table ||. Despite its simplicity. Model 
I can be used to describe and analyze the main diffu- 
sion processes: single adatom hopping, attachment, de- 
tachment and edge diffusion. The barriers obtained from 
this model can be incorporated in simulations to repro- 
duce (at least qualitatively) experimental features such 
us cluster mobility, island morphology and island density 

The model presented above describes only the gross 
features of the diffusion process. In order to get more 
quantitative results, and to better understand the impor- 
tance of the different processes, it is necessary to further 
refine the model. We will now introduce model II, in 
which the effect of NNN atoms in the initial configura- 
tion is included. The expression for the energy at the 
initial site is now 



Ei, 



^in ^ ^Enn ■ {Si + S3 + S5) 

—AEni^n ■ {So + S2 + S4 + Sq) 



(8) 



where AE^nn is the reduction of the energy due to a 
NNN bond. The energy barriers are now given by 



E'l 



Eb + AEnn ■ {S3 



-AE 



NNN 



{So 



— 5*2 — 5*6) 

S2 + S4 + So) 



(9) 



Model II accounts better for processes such as detach- 
ment, edge diffusion and vacancy diffusion, which gen- 
erally involve NNN interactions. In the distribution of 
the barriers obtained from the model, the main groups 
exhibit certain widths. Yet, they are still significantly 
narrower than the groups of the EAM barriers. This is 
due to the fact that during the hopping process adjacent 
atoms may relax within their potential well. Model II 
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accounts for these effects only on average and therefore 
gives rise to narrower groups. The values obtained as best 
fits of the model parameters to the EAM data for the dif- 
ferent metals are shown in Tabic Further refinement 
can be obtained by introducing a distinction between the 
NN bond energies at the initial four fold hollow site and 
that at the bridge site, as suggested in Ref. This 
modification, which introduces a fourth parameter into 
the model, gives only slightly better agreement with the 
EAM results. In the following Section we present a more 
effective refinement based on nonlinear interactions. 



C. Adding Non-Linear Effects 

To obtain models which provide a better fit to the 
EAM barriers, it is necessary to consider effects that are 
caused by the simultaneous interactions of the hopping 
atom with several of its neighbors. Such effects may be 
described by expressions such us SiSj or SiSjSk, which 
are equal to 1 only if all the relevant sites are occupied, 
and otherwise. Such expressions are clearly beyond 
the linear bond counting scheme of the previous Section. 
There is a large number of possible nonlinear interaction 
terms. Our analysis of the EAM calculations, however, 
indicates that two of them are most significant. The first 
term is related to the shape of the diffusion path. It 
corresponds to configurations in which sites adjacent to 
the bridge site, are occupied on both sides of the diffu- 
sion path (namely, at least one of the sites 1 and 2, as 
well as at least one of the sites 5 and 6 are occupied). 
It appears that in these cases the energy barrier is con- 
siderably higher than for configurations where sites on 
only one side of the path are occupied. This effect is due 
to the "stiffness" of the diffusion path induced by the 
attraction from two opposite directions. Even though 
there are nine different such configurations, they all con- 
tribute about the same energy difference, and hence can 
be bound to a single parameter AEopp (for opposite). 
The additional term that is now added to the expression 
for the barrier is AEopp ■ ('S'1,2 • 'S'5,6) where Sij = 1 if at 
least one of the sites i and j is occupied, and if both 
are empty. In all the metals we checked, except Cu, this 
term is much larger than the NNN bond, sometimes by 
an order of magnitude. 

The second nonlinear interaction term is smaller, and 
is comparable to the effect of NNN sites. It is related to 
the energy of the hopping atom in the initial site. The 
EAM calculations indicate that if the two nearest neigh- 
bor sites 1 and 5, that are symmetric with respect to 
the hopping direction are both occupied, then the initial 
configuration is more tightly bound. This means that 
if sites 1 and 5 are both occupied, the energy barrier is 
expected to be higher. Consequently, the corresponding 
term would be: AEg. 
tained for AE,., 



included both contributions in the same term, although 
of different physical origin, to avoid the need for a fifth 
independent parameter. 

The resulting model (model III) for the hopping energy 
barriers is 



AEnn ■ {S3 — 5*2 — 5*6) 
AEopp{Si,2 ■ S^.fi) 
AEnnn ■ {So + S4 + S1S5) 



(10) 



The values obtained from the best fit for the four pa- 
rameters E*^, AEnn, AEopp and AEnnn for the differ- 



■{Si-Sq). The fitted value ob- 
is very close to that obtained for the 
NNN binding energy E^nn, for all five metals. We thus 



ent metals are given in Table [II . There are two remarks 
to be made about Eq. (|o|). First, terms such as S1S5 
are not in contradiction to the assumption that only near- 
est and next-nearest neighbor interactions are significant. 
These terms are just a manifestation of the simultaneous 
interactions of, say the atoms in sites 5*1 and 6*5, with the 
hopping atom, of which they are both nearest-neighbors. 
Second, ^2 and 5*6 arc not included in the last term of Eq. 
(|lO|), since we found that their dominant contribution is 
in the nonlinear term. 



D. Testing the Quality of the Fit 

The quality of the fit can be viewed in Fig. |^. The 
numbering of the configurations is the decimal repre- 
sentation of the binary number n' — S3S2SeSiS5SoS4, 
where Si = 1(0) if site i is occupied (unoccupied), and 
Si is the opposite of Si. There are essentially 6 groups of 
barriers which are marked in the figures. These groups 
correspond (not necessarily in order) to the six columns 
of Table |. As can be seen, groups 11(a) and 11(b) are 
in the same energy range, and together form group II in 
Fig. H Similarly groups III (a) and Ill(b) coincide with 
group III in Fig. ^. Thus, there are actually only four 
groups as mentioned in Sec. |l|. Beyond this basic di- 
vision, there are some significant differences among the 
metals which Model III seems to handle well. The most 
important one is the effect of NNN atoms on the energy 
barrier. It can be seen from Table HI that for Ag and 
Pd this effect is almost negligible, while for Cu and Ni it 
has much greater importance. The effective NNN bind- 
ing energy for Au is even negative. Although it may be 
possible to construct models with the same number of 
parameters that would give better agreement with EAM 
results for each specific metal alone, our approach is to 
find the general characteristics of diffusion mechanisms, 
common to different substrates. The agreement between 
the EAM barriers and model III is slightly worse for Au 
than for the other metals. This may be due to substrate 
relaxation effects which are found to be more important 
in this metal, and are not accounted for in the model. 
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V. DISCUSSION AND SUMMARY 

The models presented in the previous Section help 
to identify the main physical mechanisms that deter- 
mine the activation energies of self-diffusion processes 
on FCC(OOl) metal surfaces. Although such processes 
may involve interactions with many substrate and in- 
plane atoms, they can be well described as the sum of 
few relatively simple terms. The first term is the activa- 
tion energy for hopping of an isolated atom. The main 
corrections are due to nearest-neighbor in-plane atoms at 
the initial site, as well as at the bridge site. The former 
increase the energy barrier while the latter decrease it by 
nearly the same amount. The next contribution is due 
to simultaneous presence of atoms on both sides of the 
hopping atom relative to the hopping path. This term 
is important in relatively dense environments. Its typ- 
ical value is about half that of the NN binding energy. 
A third and generally much smaller contribution consists 
of NNN bonds as well as a term associated with the si- 
multaneous presence of atoms in both sites 1 and 5 (Fig. 

Beyond the physical understanding gained by this 
analysis, the models can be used to evaluate the activa- 
tion energy of any diffusion process on the (001) surface 
of the metals discussed above. The models suggest that 
given a set of few activation energies (which can be ob- 
tained from EAM, ab-initio calculations or experiments), 
it is possible to extract the complete set of activation en- 
ergy barriers. To realize model I, for example, only two 
parameters are needed. They can be obtained e.g. from 
the activation energy for single adatom hopping, and the 
dissociation energy of a dimer. To estimate a barrier us- 
ing model II, a third parameter is needed, which is the 
NNN binding energy. This may be obtained if the mo- 
bility of a trimer is known. The fourth parameter which 
is needed for model III can be estimated from the acti- 
vation energy for detachment from an atomic step. Since 
model III provides an expression for the energy barriers, 
linear in the parameters, any four barriers that give rise 
to four linearly independent equations, are sufficient to 
determine all four parameters, and consequently all the 
other barriers. 

The possibility to construct a full set of activation en- 
ergy barriers from a relatively small set of parameters 
is especially useful for simulations. Without this knowl- 
edge, some processes have to be discarded from the sim- 
ulations as unimportant, or assigned activation energies 
which are not fully substantiated. These approaches take 
much of the power of computer simulations, and deny 
the possibility of direct quantitative confrontation with 
experimental data. Even if a list of all relevant activa- 
tion energies is available, the model can be used to check 
the self-consistency of the data. It can also help to in- 
terpret simulation results, which depend otherwise on a 
huge number of parameters. 

The models presented here apply for diffusion on flat 



surfaces and do not describe the motion up/down steps. 
Such inter-terrace moves involve a large number of possi- 
ble local enviroments, including flat steps as well as kink 
sites. We believe that the approach proposed here can 
be extended to describe these processes as well. 

In summary, we have constructed a family of models 
which describe self-diffusion on FCC (001) metal surfaces 
and tested them for Cu, Ag, Au, Ni and Pd. For each 
one of these metals, the parameters of the models were 
optimized by comparing the energy barriers to a full set of 
barriers obtained from semi-empirical potentials via the 
embedded atom method. It is found that these models, 
with at most four parameters, provide a good description 
of the hopping energy barriers on the FCC(OOl) surfaces. 

We thank G. Vidali for helpful discussions. 
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TABLE I. The hopping energy barriers for Cu, Ag, An, Ni and Pd obtained from the EAM calculations for all possible 
configurations within a 3 x 3 square around the hopping atom. The barriers are given in eV. Each number in the Table is the 
barrier for the configuration in which the occupied sites are the union of the occupied sites in the picture on top of the 
given column (indexed by ni) and on the left hand side of the given row (indexed by 72,2). Consequently, the index n specifying 
the barrier is given by n = ni + ^2 . 



TABLE IL The parameters Eq, AEnn and AEmnn of model II obtained from the best fit of the EAM barriers for Cu, 
Ag, Au, Ni and Pd. The values in parenthesis are the corresponding values for model I. R is the sum of squares defined in Eq. 
(M), for the optimized parameters. 



TABLE III. The parameters Eq, AEnn, AEnnn and AEopp of model III obtained from the best fit of the EAM barriers 
for Cu, Ag, Au, Ni and Pd. R is the sum of squares defined in Eq. (^). 
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Table II 



Metal 




.Model parameters [e\'j 




R 


Eo 


AEnn 




Cu 


0.487 (0.534) 


0.274 (0.255) 


0.027 


0.280 (0.345) 


Ag 


0.509 (0.525) 


0.204 (0.197) 


0.010 


0.450 (0.458) 


Au 


0.776 (0.752) 


0.235 (0.244) 


-0.014 


1.275 (1.292) 


Ni 


0.645 (0.697) 


0.326 (0.306) 


0.031 


0.444 (0.526) 


Pd 


0.760 (0.789) 


0.337 (0.325) 


0.017 


0.866 (0.892) 



Table III 



Metal 




Akxiel 


. parariieters[e\'j 








Eo 


AEnn 


AEnnn 
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Cu 


0.474 


0.258 


0.044 
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0.225 


0.013 


0.112 
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Au 


0.686 


0.294 


-0.017 


0.199 


0.346 


Ni 


0.610 


0.322 


0.046 


0.067 


0.141 


Pd 


0.690 


0.362 


0.028 


0.146 


0.162 
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FIG. 1. Two mechanisms for surface diffusion: (a) regular or bridge-site hopping; (b) exchange hopping. Dark spheres are 
adatoms and light spheres are substrate atoms. 

FIG. 2. Classification of all possible local environments of a hopping atom, including seven adjacent sites. Each site can be 
either occupied or unoccupied, giving rise to 2^ = 128 local environments. Sites 1, 3 and 5 are nearest neighbors of the original 
site while sites 1, 2, 5 and 6 are adjacent to the bridge site that the atom has to pass. 

FIG. 3. The distribution of activation energies barriers for Cu(OOl), Ag(OOl), Au(OOl), Ni(OOl) and Pd(OOl), obtained from 
the EAM calculations combined with the molecular statics procedure. The columns represent the number of local configurations, 
out of the 128 configurations of Fig. |^, giving rise to energy barriers in a certain energy range. Four groups of moves are identified 
in each of the five plots, and representative moves in each group are shown. 

FIG. 4. The procedure used to evaluate the binding energies. The binding energy due to two NNN bonds of Cu adatoms on 
Cu(OOl) is evaluated as the difference between the total energy Ei = —8442.043 eV of a configuration including three separate 
adatoms (a) and the total energy E2 — —8442.106 eV of a configuration including three adatoms forming two NNN bonds. 

FIG. 5. Configurations used to evaluate the NN bond energy: (a) two nearest neighbors of an adatom may be next nearest 
neighbors of each other; (b) two nearest neighbors of an atom at the bridge site, may also be nearest neighbors of each other. 



FIG. 6. Testing the additivity assumption for bond energies on the Cu(OOl) surface, (a) binding energy vs. the number of 
NNN bonds. The best linear fit yields a value of Ennn = 0.0512el/; (b) binding energy vs. number of NN bonds, where the 
binding Energy of NNN bonds is subtracted. The slope of the solid line yields -Ejvjv — 0.3246^; (c) binding energy vs. number 
of NN bonds for an atom at the bridge site, where other NN bond energies are subtracted. The best fit is Etop = 0.345el/. 
The fits include a constraint that the line passes through the origin, since an isolated adatom has zero adatom-adatom binding 
energy. 

FIG. 7. The hopping energy barriers for (a) Cu, (b) Ag, (c) Au, (d) Ni and (e) Pd, as a function of the configuration number, 
which is given by the decimal representation of the binary number n' — 5*3 ^SeSi 5*5 5*054. The solid lines are the EAM energy 
barriers, the dashed lines describe the best fits obtained for model I and the dotted lines are the best fits for model III. Model 
III is found to provide good quantitative agreement with the EAM results for most configurations. There are essentially 6 
groups of barriers in each plot, which correspond to the six columns of Table |. Groups 11(a) and 11(b) are in the same energy 
range, and together form group II in Fig. |^, while groups Ill(a) and Ill(b) coincide with group III in Fig. ^ 
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